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ABSTRACT 

We present new absolute trigonometric parallaxes and relative proper motions for nine Galactic 
Cepheid variable stars: £ Car, C Gem, (3 Dor, W Sgr, X Sgr, Y Sgr, FF Aql, T Vul, and RT Aur. We 
obtain these results with astrometric data from Fine Guidance Sensor lr, a white-light interferometer 
on Hubble Space Telescope. We find absolute parallaxes in milliseconds of arc: I Car, 2.01 ± 0.20 ; £ 
Gem, 2.78 ±0.18 ; (3 Dor, 3.14 ±0.16 ; W Sgr, 2.28 ±0.20 ; X Sgr, 3.00 ±0.18 ; Y Sgr, 2.13 ±0.29 ; FF 
Aql, 2.81 ±0.18 ; T Vul, 1.90 ±0.23 ; and RT Aur, 2.40 ±0.19 , an average ov/tt = 8%. Two stars (FF 
Aql and W Sgr) required the inclusion of binary astrometric perturbations, providing Cepheid mass 
estimates. With these parallaxes we compute absolute magnitudes in V, I, K, and Wesenheit Wyj 
bandpasses corrected for interstellar extinction and Lutz-Kclker-Hanson bias. Adding our previous 
absolute magnitude determination for S Cep, we construct Period-Luminosity relations for ten Galactic 
Cepheids. 

We compare our new Period-Luminosity relations with those adopted by several recent investiga- 
tions, including the Freedman and Sandage Ho projects. Adopting our Period-Luminosity relationship 
would tend to increase the Sandage Ho value, but leave the Freedman Ho unchanged. Comparing our 
Galactic Cepheid PLR with those derived from LMC Cepheids, we find the slopes for K and Wyj 
identical in the two galaxies within their respective errors. Our data lead to a Wyj distance mod- 
ulus for the Large Magellanic Cloud, m-M = 18.50±0.03, uncorrected for any metallicity effects. 
Applying recently derived metallicity corrections yields a corrected LMC distance modulus of (m- 
M)o=18.40±0.05. Comparing our Period-Luminosity relationship to solar-metallicity Cepheids in 
NGC 4258 results in a distance modulus, 29.28 ± 0.08, which agrees with that derived from maser 
studies. 

Subject headings: astrometry — interferometry — stars: distances — stars: individual {t Car, £ Gem, 
P Dor, W Sgr, X Sgr, Y Sgr, S Cep, FF Aql, T Vul, RT Aur) — stars: binary 
— distance scale calibration — stars: Cepheids — stars: variables — galaxies: 
individual (Large Magellanic Cloud, NGC 4258) 



1. INTRODUCTION 

Many of the methods used to determine the distances 
to remote galaxies and ultimately the size, age, and shape 
of the Universe itself depend on our knowledge of the 
distances to local objects. Among the most important 
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of these are the Cepheid variable stars. The Cepheid 
Period-Luminosity Relation (hereafter PLR) was first 
identified by Leavitt (Leavitt & Pickering 1912). This 
has led to considerable effort to determine the absolute 
magnitudes, My, of these objects, as summarized in the 
comprehensive reviews by Madore & Freedman (1992), 
Feast (1999), and Macri (2005). 

As summarized by Freedman et al. (2001) Cepheids 
are among the brightest stellar distance indicators and 
a critical initial step on the 'cosmic distance ladder'. 
These 'standard candles' are relatively young stars, found 
in abundance in spiral galaxies. For extragalactic dis- 
tance determinations many independent objects can be 
observed in a single galaxy, affording a reduction in dis- 
tance modulus error. Their large amplitudes and charac- 
teristic (sawtooth) light curve shapes facilitate their dis- 
covery and identification. Lastly, the Cepheid PLR has 
a small scatter. In the I band, the dispersion amounts to 
only -0.1 mag (Udalski et al. 1999). 

Given that the distances of all local Cepheids, ex- 
cept Polaris, are in excess of 250pc, most of the past 
absolute magnitude deter minations have used indirect 
appro ach es, for example iGroenewegen fe Oudmaiier 
(|2000f) lLanoix. Paturel. fc Gamier I d 1999ft. iFeast 
(jl997t i. lFeast fe Catchpole I (|1997h . lFeast et al. I (|1998al ) 
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and iFeast I (|1999f ). Various authors, e.g. Gieren 
et al. (1993) used Cephcid surface brightness to esti- 
mate distances and absolute magnitudes. For Cepheid 
variables, these determinations are complicated by 
dependence of the absolute magnitudes on color index 
and possibly metallicity. Only recently have relatively 
high-precision trigonometric parallaxes (— < 10%) 
been available for a very few Cepheids (the prototype, 5 
Cep and Polaris) from HIPPARCOS (jPerrvman et al I 
Il997t) . More recently we have determined the par- 
allax of 5 Cep with Fine Guidance Sensor 3 (FGS 
3) on Hub ble Space Telesc o pe (HS T) with ^ ~ 5% 
precision ([Benedict et al. I l2002b[ ). Long-baseline 
ground-based interferometry has recently provided 
radii and, through various surface brightness methods , 
distances (iNordgren et al. I l2002t iKervella et~aT[ l2004t 
lKervellall2006ll 

Our immediate goal is to determine trigonometric par- 
allaxes for an additional nine nearby fundamental mode 
Galactic Cepheid variable stars. Our target selection 
consisted in choosing the nearest Cepheids (using HIP- 
PARCOS parallaxes), covering as wide a period range 
as possible. These stars are in fact the brightest known 
Cepheids at their respective periods. Our new parallaxes 
provide distances and ultimately absolute magnitudes, 
M, in several bandpasses. Additionally, our investiga- 
tion of the astrometric reference stars provides an inde- 
pendent estimation of the line of sight extinction to each 
of these stars, a contributor to the uncertainty in the ab- 
solute magnitudes of our prime targets. These Cepheids, 
all with near solar metallicity, should be unafflicted 
by potential variations in absolute magnitude due t o 
metallicity v a riatio ns, e.g. iGroenewegen et al~l (|2004j) ; 
iMacri et al. I ([2006). Adding o ur previously de t ermine d 
absolute magnitude for S Cep iBenedict et al. I (|2002b) , 
we establish V, I, K, and Wyi Period-Luminosity Re- 
lationships (PLR) using ten Galactic Cepheids with av- 
erage metallicity, ([Fe/H])=0.02, a calibration that can 
be directly applied to external galaxies whose Cepheids 
exhibit solar metallicity. 

We describe our astrometry using one of our targets, 
£ Car, as an example throughout. This longest-period 
member of our sample provides marginal evidence for 
any possible PLR V-band non-linearity and, if included, 
anchors our PLR slopes. Additionally it is the only one 
of our sample in the period range typically used to es- 
tablish extragalactic distance moduli. Hence, its parallax 
value deserves as much external scrutiny as possible. We 
discuss (Section[2]) data acquisition and analysis; present 
the results of spectrophotometry of the astrometric ref- 
erence stars required to correct our relative parallax to 
absolute (Section [3]); derive absolute parallaxes for these 
ten Cepheid variable stars (Section 4); derive Cepheid 
absolute magnitudes (Section 5); and finally in Section 6 
we determine a number of Period-Luminosity Relations, 
briefly discuss the possibility of nonlinearity in the galac- 
tic V band PLR, discuss the distance scale ramifications 
of our results, and apply our PLR to two interesting cases 
- the LMC and NGC 4258. We summarize in Section 7. 

2. OBSERVATIONS AND DATA REDUCTION 

Nelan et al. (2003) provides an overview of the FGS in- 
strument and Benedict et al. (2002b) describe the fringe 
tracking (POS) mode astrometric capabilities of an FGS, 



along with the data acquisition and reduction strategies 
also used in the present study. We time-tag our data 
with a modified Julian Date, MJD = JD - 2400000.5, 
and abbreviate millisecond of arc, mas, throughout. 

Eleven sets of astrometric data were acquired with 
HST FGS lr for each of our nine new science targets. 
For details on our tenth, previous ly analyzed Cepheid, 
5 Cep, see IBenedict et ail (|2002bD . We obtained most 
of these eleven sets in pairs at maximum parallax fac- 
tor typically separated by a week, a strategy designed to 
protect against unanticipated HST equipment problems. 
We encountered none, obtaining 110 orbits without the 
slightest difficulty. A few single data sets were acquired 
at various minimum parallax factors to aid in separat- 
ing parallax and proper motion. Each complete data 
aggregate spans 1.49 to 1.95 years. Table Q] contains the 
epochs of observation, pulsational phase, and estimated 
B-V color index (required for the lateral color correc- 
tion discussed in Section |4~T|) for each Cepheid. The B-V 
colors are inferred from color curves constructed from 
the Cepheid photometric database 9 cited by Berdnikov 
et al. (2000). In the case of RT Aur, we supplemented 
the few data in that source with B-V values from Moffett 
k Barnes (1984), Barnes et al. (1997), and Kiss (1998). 
We adopted the periods and epochs listed by Szabados 
(1989, 1991) for these Cepheids except for I Car. Because 
the period of i Car changes unpredictably with time 
([Szabados 1 1 19891) . we derived a new period based upon 
the more recent V data in the Berdnikov et al. database. 
The Cepheids r\ Aql, £ Gem, and X Sgr also have variable 
periods, but they vary quadratically and predictably. We 
took these variations into account when computing the 
phases. 

Each individual data set required approximately 33 
minutes of spacecraft time. The data were reduced and 
calibrated as detailed in McArthur et al. (2001), Benedict 
et al. (2002a), Benedict et al. (2002b), and Soderblom 
et al. (2005). At each epoch we measured reference stars 
and the target multiple times, this to correct for intra- 
orbit drift of the type seen in the cross filter calibration 
data shown in figure 1 of Benedict et al. (2002a). A 
typical distribution of reference stars on a second gener- 
ation Digital Sky Survey R image near one of our science 
targets (£ Car) is shown in FigureQ] The somewhat elon- 
gated distribution of reference stars is forced by the shape 
of the FGS field of view and the overlap area. The orien- 
tation of each successive observation at near-maximum 
parallax factor changes by 180°, mandated by HST solar 
panel illumination constraints. 

Data are downloaded from the HST archive and passed 
through a pipeline processing system. This pipeline ex- 
tracts the astrometry measurements (typically one to two 
minutes of fringe position information acquired at a 40 
Hz rate, which yields several thousand discrete measure- 
ments), extracts the median (which we have found to be 
the optimum estimator), corrects for the Optical Field 
Angle Distortion (McArthur et al. 2002), and attaches 
all required time tags and parallax factors. 

Table [5] collects measured properties for our tar- 
get Cepheids, including pulsational period, log of 
that period, (V) ,(I), (K) (B-V), E(B-V ) Ay, and 
Ak ■ Photometry is from IGroenewegen I (|1999f ) and 

9 |http : //ftp . sai .msu . su /groups/cluster/CEP/PHE 
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iBerdnikov et al. I (|1996f ). The (K) values for 6 Cep 
and T Vul were corrected following Berdnikov (2006, 
private communication). Cepheid (K) is in the CIT 
system. The (I) is in the Cousins system. All 
reddening values are either derived from our refer- 
ence star photometry or adopted from those listed 
in the David Dunlap Observatory C epheid database 
(|Fernie. Evans. Beattie. fc Seager Ill995[ ) . Our reddening 
selection criterion is discussed in Section 5. 

3. SPECTROPHOTOMETRIC PARALLAXES OF THE 
ASTROMETRIC REFERENCE STARS 

The following review of our astrometric and spec- 
trophotometric techniques uses the £ Car field as an 
example. Given that £ Car has the longest period in our 
sample, it may have a significant effect on the slopes of 
the PLR we eventually construct. It also has a period 
most like that of extragalactic Cepheids used in distance 
determination. Because the parallaxes determined for 
the Cepheids will be measured with respect to reference 
frame stars which have their own parallaxes, we must ei- 
ther apply a statistically derived correction from relative 
to absolute parallax (Van Altena, Lee & Hofleit 1995, 
hereafter YPC95) or estimate the absolute parallaxes of 
the reference frame stars listed in Table [5] In principle, 
the colors, spectral type, and luminosity class of a star 
can be used to estimate the absolute magnitude, My, 
and V-band absorption, Ay. The absolute parallax is 
then simply, 



The luminosity class is generally more difficult to esti- 
mate than the spectral type (temperature class). How- 
ever, the derived absolute magnitudes are critically de- 
pendent on the luminosity class. As a consequence we 
use as much additional information as possible in an at- 
tempt to confirm the luminosity classes. Specifically, we 
obtain 2MASS 1 photometry and UCAC2 proper motions 
(Zacharias et al. 2004) for a one degree square field cen- 
tered on each science target, and iteratively employ the 
technique of reduced proper motion (Yong & Lambert 
2003, Gould & Morgan 2003, Ciardi 2004) to confirm 
our giant/dwarf classifications. 

3.1. Reference Star Photometry 

Our band passes for reference star photometry in- 
clude: BVI (from recent measurements with the New 
Mexico State University lm telescope for the northern 
Cepheids, and from the South African Astronomical Ob- 
servatory (SAAO) lm for the southern Cepheids) and 
JHK (from 2MASS). For reference star spectrophotomet- 
ric parallaxes only, the 2MASS JHK have been trans- 
formed to the Bessell & Brett (1988) system using the 
transformations provided in Carpenter (2001). Table [3] 
lists BVIJHK photometry for targets and reference stars 
bright enough to have 2MASS meas urements. In addi- 
tion Washington-DDO photo metry (jPaltoglou k, Bell~l 
[1991 IMaiewski et al I l200l was used to confirm the 
luminosity classifications for the later spectral type ref- 
erence stars. 

The Two Micron All Sky Survey is a joint project of the Uni- 
versity of Massachusetts and the Infrared Processing and Analysis 
Center / California Institute of Technology 



3.2. Reference Star Spectroscopy 

The spectra from which we estimated I Car reference 
star spectral type and luminosity class come from the 
South African Astronomical Observatory (SAAO) 1.9m 
telescope. Spectral classifications for the [3 Dor and X, 
Y, and W Sgr fields were also provided by SAAO. The 
SAAO resolution was 3.5 A/ (FWHM) with wavelength 
coverage from 3750 A< A < 5500 A. Spectroscopic clas- 
sification of the reference stars in the fields of RT Aur 
and Q Gem was accomplished using data obtained with 
the Double Imaging Spectrograph (DIS) on the Apache 
Point Observatory 3.5 m telescope . We used the high 
resolution gratings, delivering a dispersion of 0.62 A/pix, 
and covering the wavelength range of 3864 < A < 5158 A. 
Spectroscopy of the reference stars in the fields of Y Sgr, 
FF Aql, ij Aql, and T Vul was obtained using the R-C 
Spectrograph on the KPNO 4 m. The "t2kb" detector 
with grating "#47" was used to deliver a dispersion of 
0.72 A/pix, covering the wavelength range 3633 < A < 
5713 A. Classifications used a combination of template 
matching and line ratios. Spectral types for the stars are 
generally better than ±2 subclasses. 

3.3. Interstellar Extinction 

To determine interstellar extinction we first plot the 
reference stars on a J-K vs. V-K color-color diagram. 
A comparison of the relationships between spectral type 
and intrinsic color against those we measured provides 
an estimate of reddening. Figure [2] contains the £ Car J- 
K vs V-K color-color diagram and reddening vector for 
Ay = 1.0. Also plotted are mappings between spec- 
tral type and luminosity class V and III from Bessell 
& Brett (1988) and Cox (2000). Figure^ along with the 
estimated spectral types, provides an indication of the 
reddening for each reference star. 

Assuming an R = 3.1 Galactic reddening law (Savage 
& Mathis 1979), we derive Ay values by comparing the 
measured colors (Table [3] ) with intrinsic (V-K)o, (B- 
V) , (U-B) , (J-K)o, and (V-I) , colors from Bessell & 
Brett (1988) and Cox (2000). We estimate Ay from Ay 
= l.lE(V-K) = 5.8E(J-K) = 2.77E(U-B) = 3.1E(B-V) 
= 2.26E(V-I), where the ratios of total to selective ex- 
tinction were derived from the Savage & Mathis (1979) 
reddening law and a reddening estim ate in the direction 
of £ Car from ISchlegel et alj <|l998), via NED 3 . All re- 
sulting Ay are collected in Table |4] We then calculate a 
field wide average Ay to be used in equation 1. For the 
£ Car field (Ay) = 0.52±0.06 magnitude. In this case 
our independent determination is in good agreement with 
the David Dunlap Observatory online Galactic Cepheid 
database 4 , which averages seven measurements of color 
excess to obtain (E(B-V)) = 0.163±0.017 , or (A v ) = 
0.51±0.05. 

Using the £ Car field as an example, we find that the 
technique of reduced proper motions can provide a pos- 
sible confirmation of reference star estimated luminosity 
classes. The precision of existing proper motions for all 
the reference stars was ~5 mas y _1 , only suggesting dis- 
crimination between giants and dwarfs. Typical errors 

2 The Apache Point Observatory 3.5 m telescope is owned and 
operated by the Astrophysical Research Consortium. 

3 NASA/IPAC Extragalactic Database 

4 http: / /www. astro. utoronto.ca/DDO/research/cepheids/cephcids. html 
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on Hr, a parameter equivalent to absolute magnitude, 
My, were about a magnitude. Nonetheless, a reduced 
proper motion diagram did suggest that ref-4, -5, and -8 
are not dwarf stars. They are considerably redder in J-K 
than the other stars in the present program classified as 
dwarfs. Giants are typically redder in J-K than dwarfs 
for a given spectral type (Cox 2000). Our luminosity 
class uncertainty is reflected in the input spectrophoto- 
metric parallax errors (Table [5]). We will revisit this ad- 
ditional test in Section ^. f 1 once we have higher precision 
proper motions obtained from our modeling. 

3.4. Estimated Reference Frame Absolute Parallaxes 

We derive absolute parallaxes f or each reference star 
using My values from lCoxH (|2000h and the (Ay) derived 
from the photometry. Our adopted errors for (m-M)o are 
0.5 mag for all reference stars. This error includes un- 
certainties in (Ay) and the spectral types used to esti- 
mate My. Our reference star parallax estimations from 
Equation f are listed in Table [5l For the £ Car field 
individually, no reference star absolute parallax is bet- 
ter determined than — = 23%. The average absolute 
parallax for the reference frame is (ir a b s ) = 0.85 mas. 
We compare this to the correction to absolute parallax 
discussed and presented in YPC95 (section 3.2, fig. 2). 
Entering YPC95, fig. 2, with the I Car Galactic latitude, 
1 = -7°, and average magnitude for the reference frame, 
(V r ef)= 13.0, we obtain a correction to absolute of 1 
mas. This gives us confidence in our spectrophotometric 
determination of the correction to absolute parallax. As 
in past investigations we prefer to introduce into our re- 
duction model our spectrophotmetrically estimated ref- 
erence star parallaxes as observations with error. The 
use of spectrophotometric parallaxes offers a more direct 
(less Galaxy model-dependent) way of determining the 
reference star absolute parallaxes. 

4. ABSOLUTE PARALLAXES OF GALACTIC CEPHEIDS 
4.1. The Astrometric Model 

With the positions measured by FGS lr we determine 
the scale, rotation, and offset "plate constants" relative 
to an arbitrarily adopted constraint epoch (the so-called 
"master plate" ) for each observation set (the multiple ob- 
servtions of reference stars and Cepheid targets acquired 
at each epoch listed in Table 1). The rotation to the sky 
of the master plate is initially set at a value provided by 
the HST ground system. The mJD of each observation 
set is listed in Table [TJ along with a Cepheid B-V esti- 
mated from a phased light curve. Our £ Car reference 
frame contains 6 stars. All the Cepheid primary science 
targets, including £ Car, are bright enough to require the 
use of the FGS neutral density filter. Hence, we use the 
modeling approach outlined in Benedict et al. (2002b), 
with corrections for both cross-filter and lateral color po- 
sitional shifts, using values specific to FGS lr determined 
from previous calibration observations with that FGS. 

We employ GaussFit (Jefferys et al. 1988) to mini- 
mize x 2 : our model goodncss-of-fit metric. GaussFit 
has a number of features, including a complete program- 
ming language designed especially to formulate estima- 
tion problems, a built-in compiler and interpreter to sup- 
port the programming language, and a built-in algebraic 
manipulator for calculating the required partial deriva- 



tives analytically. The program and sample models are 
freely available 5 . 

The solved equations of condition for the I Car field 
are: 

x' = x + lc x (B ~ V)~AXFx (2) 

y' = y + lCy(B- V)-AXFy (3) 

£ = Ax' + By' + C - HxAt - P a ir x (4) 

r] = Dx + Ey' + F - (iyAt - P S TT y (5) 

where x and y are the measured coordinates from HST: 
lc x and lc y are the lateral color corrections; AXFx and 
AXFy are the cross filter corrections in x and y, applied 
only to the observations of each Cepheid; and B — V are 
the B-V colors of each star. A, B, D, and E are scale 
and rotation plate constants, C and F are offsets; \x x and 
\x v are proper motions; At is the epoch difference from 
the mean epoch; P a and P$ are parallax factors; and n x 
and TT y are the parallaxes in x and y. x' and y' are FGS 
positions corrected for lateral color and cross-filter shifts. 
£ and 77 are relative positions in arcseconds. We obtain 
the parallax factors from a JPL Earth orbit predictor 
(Standish 1990), upgraded to version DE405. 

There are additional equations of condition relating an 
initial value (an observation with associated error) and 
final parameter value. There is one such equation in the 
model for each parameter of interest: reference star and 
target color index, proper motion, and (excepting the 
Cepheid target) spectrophotometric parallax. Through 
these additional equations of condition the y 2 minimiza- 
tion process is allowed to adjust parameter values by 
amounts constrained by the input errors. We also sim- 
ilarly adjust the lateral color parameters, master plate 
roll, and cross filter parameters. The end results are the 
final values of the parameters of intertest. In this quasi- 
Bayesian approach prior knowledge is input as an obser- 
vation with associated error, not as a hardwired quantity 
known to infinite precision. 

For example input proper motion values have typical 
errors of 4-6 mas y _1 for each coordinate. Final proper 
motion values and errors obtained from our modeling 
of HST data for the £ Car field are listed in Table 
Adjustments to the proper motion estimates required to 
minimize x 2 averaged 3 mas yr _1 . For completeness, 
transverse velocities, given our final parallaxes, are listed 
in Table | As a final test of the quality of our prior 
knowledge of reference star luminosity class listed in Ta- 
ble [5l we employ the technique of reduced proper mo- 
tions. We obtain proper motion and J, K photometry 
from UCAC2 and 2MASS for a |°x |° field centered 
on £ Car. Figure [3] shows H^- = K + 51og(/x) plotted 
against J-K color index for 436 stars. If all stars had the 
same transverse velocities, Figure [3] would be equivalent 
to an HR diagram. £ Car and reference stars are plotted 
as ID numbers from Table [JJ I Car is T' in Figured 
With our precise proper motions (Table [7J errors in H^- 
are now ~ 0.3 magnitude. Reference stars ref-4, -5, and 
-8 remain clearly separated from the others, supporting 
their classification as giants. 

We stress that for no Cepheid in our program was a 
previously measured parallax used as prior knowledge 

5 http:/ /clyde. as. utexas.edu/Software. html 
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and entered as an observation with error. Only reference 
star prior knowledge was so employed. Our Cepheid par- 
allax results are blind to previous parallax measures from 
HIPPARCOS and/or parallaxes from surface brightness 
estimates. 

4.2. Assessing Reference Frame Residuals 

The Optical Field Angle Distortion calibration 
(McArthur et al. 2002) reduces as-built HST telescope 
and FGS lr distortions with amplitude ~ 1" to below 
2 mas over much of the FGS lr field of regard. From 
histograms of the i Car field astrometric residuals (Fig- 
ure U) we conclude that we have obtained satisfactory 
correction. The resulting reference frame 'catalog' in £ 
and rj standard coordinates (Table [6]) was determined 
with average position errors (o^) = 0.50 and (ov,) = 0.62 
mas. 

To determine if there might be unmodeled - but pos- 
sibly correctable - systematic effects at the 1 mas level, 
we plotted reference frame X and Y residuals against 
a number of spacecraft, instrumental, and astronomical 
parameters. These included X, Y position within our to- 
tal field of view; radial distance from the field of view 
center; reference star V magnitude and B-V color; and 
epoch of observation. We saw no obvious trends. 

4.3. Absolute Parallaxes of the Cepheids 

For the I Car and Y Sgr fields we reduced the num- 
ber of modeling coefficients in equations 3 and 4 to 
four, as done for our previous work on the Pleiades 
(|Soderblom et ai"1l2005f ). We constrained the fit to have 
a single scale term by imposing D = -B and E = A. 
Final model selection for all fields was based on refer- 
ence star placement relative to the target, total number 
of reference stars, reduced x 2 (x 2 /DOF, where DOF = 
degrees of freedom), and parallax error. Absolute paral- 
laxes, relative proper motions, and transverse velocities 
for I Car and associated reference stars are collected in 
Tables [7] and [H Parallaxes for all Cepheids are collected 
in Table HU 

All our Cepheid parallaxes directly rely on the esti- 
mates of reference star parallaxes. Should anyone wish 
to independently verify our results, the reference stars 
used in this study are all identified in archival material 6 
held at the Space Telescope Science Institute. Adopted 
reference star spectral types and the parallaxes resulting 
from our modeling are listed in Table [9] Similar data for 
the S C ep reference stars can be found in lBenedict et all 
(|2002bH . 

4.4. HST Parallax Accuracy 

Our parallax precision, an indication of our internal, 
random error, is ~ 0.2 mas. To assess our accuracy, or 
external error, we have compared (Benedict et al. 2002b, 
Soderblom et al. 2005) our parallaxes with results from 
independent measurements from HIPPARCOS (Perry- 
man et al. 1997). Other than for the Pleiades (Soderblom 
et al. 2005) , we have no large systematic differences with 
HIPPARCOS for any objects with ^ <10%. The next 
significant improvement in geometrical parallaxes for 
Cepheids will come fro m the space-ba sed, all-sky astrom- 
etry missions GAIA (jMignardl [2005T ) and SIM (jUnwinl 

6 [EJtp: / / www.stsci.edu/observing/ phase2-public / 9879. pro 



I2005T) with - 10 /larcsec precision parallaxes. Final re- 
sults are expected by the end of the next decade. 

4.5. The Binary Cepheids 

Many of our target Cepheids have companions discov- 
ered spectroscopically with IUE (c. /. Eva ns 1995). 
Two of these, W Sgr (|Petterson et al. I l2004h and FF 
Aql ([Evans et al.lll990l ). have published spectroscopic or- 
bital elements. For these two targets we introduced the 
known orbital elements as observations with error and 
solve for inclina t ion and perturbation size as outlined in 
iBenedict et al. I (|2002b), using equations 6 and 7 from 
that paper. 

Our results for W Sgr and FF Aql are summarized in 
Table flOl With the perturbation orbit semimajor axis, a, 
the measured inclination, i, and an estimate of the sec- 
ondary mass, we can estimate the mass of each Cepheid. 
The secondary mass is estimated from the spectral typ e 
and a recent Mass-Luminosity relationship (|Henrvl l2004') . 
An improvement in the W Sgr mass is expected shortly, 
once the secondary spectral type is more tightly con- 
strained (Evans & Massa 2007). The major contributor 
to the FF Aql mass error is the parallax error. 

One of our original targets, r\ Aql, is thought to be a 
binary from IUE spectra (B9.8 companion, Evans 1991). 
As shown in Table fTOl we have successfully included per- 
turbation orbits for two other Cepheids, FF Aql and 
W Sgr, simultaneously solving for parallax, proper mo- 
tion, inclination, and perturbation orbit semimajor axis. 
However, spectroscopic orbital parameters are fairly well 
known for those stars, definitely not the case for r\ Aql. 
With no period, eccentricity, or periastron timing con- 
straints from previous radial velocity observations, and 
effectively only five distinct epochs of astrometry, we 
cannot determine a perturbation orbit. Our astrome- 
try is clearly affected by an, as yet unmodelable, motion. 
Therefore, r\ Aql cannot be included in this analysis. Ul- 
timately, additional HST observations may serve to char- 
acterize the (likely face-on) perturbation orbit, resulting 
in a usable parallax. 

5. THE ABSOLUTE MAGNITUDES OF THE CEPHEIDS 

When using a trigonometric parallax to estimate the 
absolute magnitude of a star, a correction should be 
made for the Lutz-Kelker bias (|Lutz fc Kelker|[l973f ) as 
modified by Hanson (1979). We justify the application 
of Lutz-Kelker- Hanson (LKH ) with an appeal to Bayes 
Theorem. See iBarnes et al. I (|2003f ). section 4, for an 
accessible introduction to Bayes Theorem as applied to 
astronomy. Invoking Bayes Theorem to assist with gener- 
ating absolute magnitudes from our Cepheid parallaxes, 
one would say, "what is the probability that a star from 
this population with this position would have parallax 7r 
(as a function of ir) , given that we haven't yet measured 
7r?" In practice one would use the space distribution of 
the population to which the star presumably belongs. 
This space distribution is built into the prior p(7r) for 7r, 
and used to determine 

p(n\ir observe d&,K) ~ p{-K observed \TT&K)p(ii\K) (6) 

where K is prior knowledge about the space distribution 
of the class of stars in question and '&' is an 'and' opera- 
tor. The function observed^ & K) is the standard like- 
lihood function, usually a gaussian normal with variance 



6 



Benedict et al. 



ov The "standard" L-K correction has p(7r| K) ~ 7r~ 4 . 
Looking at a star in a disk population close to the galac- 
tic plane requires 7r~ 3 (ignoring spiral structure), which 
is the prior we use. The LKH bias is proportional to 
(ov/7r) 2 . Presuming that all Cepheids in Table 2 be- 
long to the same class of object (evolved Main Sequence 
stars), we scale the LKH correction determined in Bene- 
dict et al. (2002b) for 8 Cep and obtain the LKH bias 
corrections listed in Table [HJ For I Car we find LKH = 
-0.08 magnitude. The average LKH bias correction for all 
Cepheids in this study was -0.06 magnitude. We identify 
the choice of prior for this bias correction as a possible 
contributor to systematic errors in the zero points of our 
PLR, at the 0.01 magnitude level. 

With (V)= 3.724 (Tabled and given the absolute par- 
allax, 2.01 ± 0.20 mas from Section 14. 3[ we determine 
a distance modulus for t Car. From Table 2] (Section 
!3-3[) we obtain a derived field-average absorption, (Ay) 
= 0.52. With this (Ay), the measured distance to I Car, 
and the LKH correction we obtain My = —5.35 ± 0.22 
and a corrected true distance modulus, (m-M)o = 8.56. 
The My error has increased slightly by combining the 
(Ay) error and the raw distance modulus error in quadra- 
ture. The Mk values in Table QTJ have slightly lower er- 
rors because the (A^) values are lower with correspond- 
ingly lower errors to add in quadrature. The Wesenheit 
magnitude, Wyj , listed in Table fTTI is the prescription of 
iFreedman et~aT\ if^OOlh . W V i = V '- E (v V _ r) (V-I). For 
the reddening law adopted by them Wvi = V — 2A5(V — 
I)- 

Results, including all proper motions, and absorption- 
and LKH bias-corrected absolute magnitudes, for the 
Cepheids in our program (except rj Aql) are collected 
in Table QT] In half the cases the reddening values we 
derived from our reference star photometry agreed with 
that listed in the David Dunlap Observatory Cepheid 
database. For 5 Cep and X Sgr we adopted a reddening 
derived as described in Section 3.3 because the photome- 
try showed very consistent star-to-star reddening. For (3 
Dor, Y Sgr, and FF Aql we adopted the DDO color ex- 
cess and an absorption Ay=3.1E(B-V) because the star- 
to-star reddening indicated extremely patchy absorption. 
Adopted absorption in the V and K band are listed in 
Table H 

6. PERIOD-LUMINOSITY RELATIONS, DISTANCE SCALE 
IMPLICATIONS, AND APPLICATIONS 

6.1. Period- Luminosity Relations from HST Parallaxes 

Plotting the absorption and LKH bias corrected V, I, 
K and Wesenheit absolute magnitudes, My, M], Mk, 
and M\y(yj), from Table ITTI against the logarithm of the 
period (logP, Table [2|) we obtain the Period-Luminosity 
relationships (PLR) contained in Figure [5] We param- 
eterize all PLR as Mx = a + b(logP — 1). Hence, the 
zero-points are for a Cepheid with logP = 1. Our inter- 
cepts and slopes (B07) with la errors are collected in Ta- 
ble!^ along with other recent determinations; Freedman 
et al (2001, F01), Sandage et al (2004, S04), and Barnes 
et al. (2003, B03). Note that th e Sandage M^yj) wa s 
derived from their V and I PLR (ISandage et al I (|2004h . 
cqs. 17 and 18), using W VI = V-2.52(V-I). Adopting 
the Freedman reddening coefficient (2.45) would change 
the slope of the Sandage Wyi PLR only by +0.02 mag- 



nitude. 

The standard deviation of our residuals (Figure [5]) are 
0.10 magnitude for My and Mj, and 0.09 magnitude for 
Mk and M W ( VI y In each case the largest residual is 
that of W Sgr. Note that the determination of the W 
Sgr parallax was complicated by the inclusion of a binary 
perturbation orbit. However, excluding W Sgr from the 
fit of, for example, M^ , changes the slope and intercept 
by less than 0.01 magnitude. 

Given the diversity of opinion regarding the applica- 
bility of LKH bias corrections (e.g. Smith 2003), even 
among the present author list, one of us (MF) suggested 
the following. The absolute magnitude error depends on 
the fractional parallax error, o-^/n. When forming the 
PLR in Figure [5] a star, which has by chance an over- 
estimated parallax, will have greater weight in the solu- 
tion than the same star with an underestimated parallax. 
Feast (1998b), table 1, presents an extreme example of 
this. A correction based on a^/n seems required. We 
first fit a PLR with absolute mags uncorrected for LKH 
bias and weight the stars by (^) 2 . From the deviations 
of each star from this PLR we deduce the parallax (iri) 
which each star would have to have to fall on the PLR. 
We then redo the PLR, weighting the uncorrected ab- 
solute mags by (■^) 2 ■ Changes in the PLR slope and 
intercept are less than 0.005 within three iterations. Ta- 
ble lists PLR slopes and intercepts obtained with this 
particular weighting scheme as B07f. We note that the 
B07f slopes and intercepts agree within their respective 
errors with those obtained employing LKH bias correc- 
tions (B07). 

6.2. Pulsation Modes and Shocks 

There is always a possibility that a Cepheid will be 
pulsating in an overtone, especially at shorter periods. 
The ratio of the fundamental period Pq to the overtone 
period Pi is given by : 

Pi/P = 0.720 - 0.0271ogP (7) 

(e.g. Alcock et al. 1995). Thus for the My - logP 
plot found in Figure [5] overtones will lie about 0.4 mag- 
nitude above fundamental pulsators at a given period. 
As Figure [5] shows there is no evidence that any of the 
Cepheids in our sample deviate by this amount from a 
single relation. This is particularly interesting in the case 
of FF Aql. This low amplitude (AV = 0.33) Cepheid 
has been classed as an overtone pulsator on the basis 
of fourier analyses of both the light curve (Antonello 
et al 1990) and the velocity curve (Kienzle et al 1999). 
It was nevertheless classed as a fundamental pulsator by 
Sachkov (1997) from a Baade-Wesselink type radius es- 
timate. None of the other Cepheids in our sample have 
been suggested to be overtone pulsators so far as we are 
aware, although ^ Gem has a relatively low amplitude 
(AV = 0.49). As shown in Table fTTI our absolute magni- 
tude for FF Aql (logP = 0.6504) is My = -3.05 ± 0.15. 
This is 0.15 magnitude fainter than the mean relation 
(Figure [5|) . If it were an overtone pulsator it would be 
0.53 magnitude fainter than expected. We therefore con- 
clude that it is not an overtone pulsator. This results 
suggests that despite the rather clear division of light and 
velocity fourier coefficients of Cepheids into two group- 
ings, this may not always correspond to a division into 
fundamental and overtone pulsators. 
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Mathias et al. (2006) find that X Sgr is apparently 
unusual in showing evidence of multiple shock waves in 
its atmosphere. It is therefore worth noting that this 
Cepheid appears quite normal as far as its position in 
any of our PLR is concerned. 

6.3. Distance Scale Implications 

In this section we use two methods to compare our 
new parallaxes and our derived PLR with those from 
previous investigations. The first approach is simply to 
compare zero-points and slopes of various PLR. These 
will be flagged as DC (direct com parison). Th e second 
approach uses reduced parallaxes (jFeast II2002D to solve 
for independent zero-points, which are then compared 
with our new values listed in Table Q21 This approach 
is denoted RP (reduced parallaxes). In the following 
we carry out these comparisons with two populations of 
Cepheids, first Galactic Cepheids, then Cepheids in the 
Large Magellanic Cloud (LMC). 

6.3.1. Comparison with Other Galactic PLR 

6.3.1.1 - Barnes et al. (2003): develop and describe a 
Bayesian statistical analysis to solve the surface bright- 
ness equations for Cepheid distances and stellar proper- 
ties. Their analysis averages over the probabilities as- 
sociated with several models rather than attempting to 
pick the best model from several possible models. They 
obtain a PLR using a sample of 13 Galactic Cepheids. 

Before comparing their V band PLR with ours, we can 
compare our parallax for T Vul with that determined 
by Barnes et al. (2003) using a Bayesian solution in 
the visual surface brightness technique. As discussed by 
Barnes et al. the surface brightness technique determines 
a quasi-geometric parallax rather than a distance. Their 
table 7 gives a parallax for T Vul of 1.65 ± 0.11 mas, 
compared to our HST parallax of 1.90 ± 0.23 mas. 

However, the surface brightness parallax contains an 
unknown systematic uncertainty, because it depends 
upon the adopted factor for conversion from radial ve- 
locity to pulsational velocity, normally denoted p. (The 
larger the value of p, the smaller the parallax computed.) 
The value adopted by Barnes et al. for T Vul is p — 1.37. 
Recent values in the literature, appropriate for periods 
near that of T Vul, range from p = 1.27±0.06 (Merand et 
al. 2005) top = 1.48 ±0.21 (Gieren et al. 2005). We can 
use the HST parallax to infer a quasi-geometrical value 
of p by demanding that the surface brightness parallax 
match the HST parallax. The result is p = 1.19 ± 0.16. 
This value is consistent with the geometrically deter- 
mined value of Merand et al., with the canonical value 
p = 1.36 that is often used for Cepheids, and with the 
value p = 1.48 from Gieren et al. (2005) given the mu- 
tual uncertainties. Even though the uncertainty is larger 
than we would like, our value is only the second geomet- 
rically or quasi-geometrically determined value of p after 
that of Merand et al. (2005). 

DC A comparison of the Barnes et al. PLR in V (based 
on thirteen Galactic Cepheids) with that for our ten- 
Cepheid solution is quite satisfactory. As shown in Ta- 
ble [T3] the zero points and the slopes both agree within 
1(7. The agreement would have been slightly better had 
the reddening law chosen by Barnes et al. been the same 
as adopted in this work. Their choice of R ^3.35 leads 
to larger values of A v than does the law adopted here. 



Adjusting the surface brightness My values to our red- 
dening law would change the slope determined by Barnes 
et al. about 0.5c closer to the slope determined in the 
present work. 

6.3.1.2 - Freedman and Sandage: An ultimate goal 
of many workers has been to use Cepheids to establish 
the Hubble Constant, Hq. Two major groups are re- 
cently or presently involved in this effort and their results 
are summarized in Freedman et al. (2001) and Sandage 
et al. (2006). Here we investigate how our new parallaxes 
test some of their basic assumptions. 

Both groups effectively use 'reddening free' (Wvi type) 
relations. However they differ, among other things, by 
use of different reddening laws and therefore different 
color coefficients in their relations. The relation used by 
Freedman et al. (2001) is equation 5 of their section 3.3 
and is derived from the OGLE LMC work on Cepheids 
(Udalski et al. 1999). With the other equations in that 
section it can be written 

Mod = V + 3.255(logP - 1) - 2.45(^ — I) + F (8) 

where the Freedman zero-point, F = +5.899+A, is based 
on an LMC distance modulus of 18.50 and a correction, 
A, for any metallicity difference with the LMC. The re- 
lation between A and metallicity adopted by Freedman 
et al. leads to A = +0.08 in the case of metal normal 
(Galactic) Cepheids. 

Sandage et al. adopt different relations for the LMC 
and our Galaxy. For our Galaxy they adopt 

Mod = V + 3.746(Zo 3 P - 1) - 2.523(V - I) + S (9) 

where the Sandage zero-point S = +5.959 is based 
on data from Baade-Wesselink type analyses and from 
Cepheids in Galactic clusters. 

RP Our first test uses the method of reduced paral- 
laxes outlined in Feast (2002, eqs 1, 2, and 3) to esti- 
mate the zero-points, F and S, in equations 8 and 9. 
The results are shown in Table [T3l together with the val- 
ues adopted by the two groups. In our reductions we 
have assumed that all the uncertainty is in the paral- 
laxes. The small scatter about the PLR relations dis- 
cussed in this paper suggest that other sources of uncer- 
tainty are small. This assumes (equation 3 of Feast 2002) 
that the uncertainty in the magnitude and the intrinsic 
scatter in the adopted relation are small. If they were 
significant they would decrease the value of the derived 
constants in the above Equations 9 and 10 (decrease the 
absolute brightness of the Cepheids) due to a change in 
the relative weights of the stars. For instance in the 
case of S this changes from +5.96 to +5.93 if each of 
the above uncertainties was 0.07 mag and to +5.91 if 
they were both 0.1 mag. The results in Table [TBI show 
that the reduced parallax calculation gives a zero-point 
for the Sandage Galactic relation in agreement with that 
adopted by them, whereas the Freedman zero-point is 
0.16 mag smaller than the one used by them. 

DC For our second test we calculate parallaxes for our 
10 Cepheids using Equations 8 and 9 with the Freedman 
and Sandage zero-points. These 'post-diction' parallaxes 
are listed in Table [T4l The results in Table [T4l confirm the 
differences found in the RP test. The small differences 
between Tables [T3] and Q3] are due to the fact that the 
DC test (Table Q3]) was carried out using unweighted 
quantities 
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In summary, using RP and our Table [TT] parallaxes 
to derive the zero-points of Equations 8 and 9, we find 
an agreement with Sandage et al. (Eq 9) within the er- 
rors but a difference of 0.16 ± 0.04 from Freedman et al. 
(Eq 8). However most HST-based work on extragalactic 
Cepheids has been heavily weighted to the longer periods, 
true both of the Freedman et al. and Sandage et al. pro- 
grams. Thus it is important to note (Table [M]) that for 
our longest period Cepheid, I Car, the Freedman et al. re- 
lation predicts a parallax in better agreement with ours 
than does that of Sandage et al. . This DC result is 
shown in a slightly different way in Fig [5] where we plot 
both the Freedman et al. and the Sandage et al. Wesen- 
heit M W ( V i) PLR together with those derived from our 
data. We allow both the coefficient of (logP - 1) and the 
zero-point to vary, but adopted the relevant color coeffi- 
cient. These plots show that at the longer periods (log P 
> f ) relevant to much extragalactic work the Freedman 
et al. relation lies close to our best estimate and im- 
plies little change to their derived Ho, despite its being 
inconsistent with our data at shorter periods. On the 
other hand using our our M w ^ V i) PLR (Fig O would 
tend to increase the Sandage estimate of Ho, at least 
where it depends on galaxies of near solar metallicity. 
This result depends crucially on £ Car. It would clearly 
be important to strengthen the long period calibration 
by obtaining additional high-precision parallaxes of long- 
period Cepheids. 

6.3.2. Comparison with LMC Cepheid PLR 

To carry out direct comparisons of Galactic and LMC 
PLR (DC) it is necessary to compare PLR with equal 
slopes, because the logP=l intercept depends on PLR 
slope. Where necessary we have refit our PLR, constrain- 
ing the slope to those established for LMC Cepheids, to 
re-determine zero-points, all uncorrected for metallicity 
effects. 

6.3.2.1 - LMC V-ban d PL R Non-linearity: 

INgeow et all (|2005l h iNgeowl (|2006l ); iNgeow k Kanburl 
(2006) offer further evidence for a possible change of 
slope of the V-band PLR at logP^l for the Cepheids 
in the LMC first noted by Sandage et al. . However, 
they find no such slope change in any reddening-free 
Wesenheit magnitude. DC With our small scatter 
V-band PLR do Galactic Cepheids exhibit a similar 
non-linearity? Fitting only those seven Cepheids in 
Figure [5] with logP < 1 for the V-band PLR results in 
a slope with a significant error. Our fit to that period- 
restricted subset is shown in Figure [71 and demonstrates 
that our sample, containing only one Cepheid with 
logP well in excess of unity, is too small to offer solid 
evidence for a V-band PLR slope change similar to that 
found in the LMC by Ngeow et al. . We note that I 
Car lies only ~ 2cr below the relationship line. We also 
note that we can obtain a V-band PLR slope only la 
different from the full sample by retaining £ Gem, and 
(3 Dor in the sample. Suspected non-linearity in V rests 
entirely on I Car. Because the PLR in V is a collapsed 
period-luminosity-color relation, it has a finite width in 
V (e.g., Caldwell k Coulson 1986). The location of I 
Car below the PLR fitted to shorter period stars could 
be a result of this finite width. 

6.3.2.2 - Gieren et al. (2005): contains absolute mag- 
nitude information for thirteen LMC Cepheids derived 



using surface brightness methods. In addition to high- 
precision photometry, ten of the thirtee n have metal- 
licity measures (([Fe/H])=-0.46). Because [Gieren et all 
( 2005) contained none of the observed intensity-averaged 
apparent magnitudes, J. Storm, a co-author on that pa- 
per, kindly supplied their (V), (I), and (K) for these 
selected LMC Cepheids. The (V), (I) and (K) were not 
corrected for LMC tilt, but are corrected for absorption 
(A v = 3.1E(B-V), A K = 34E(B-V)), using the E(B-V) 
from iGieren et al. I (|2005h ). The V and K PLR, along 
with a Wesenheit W VI PLR (W V/ =V-2.45(V-I)), are 
shown in Figure [5] 

DC Comparing the Gieren et al. LMC PLR Figure [5] 
PLR with our Galactic PLR in Figure [5] we note satisfac- 
tory agreement in the slopes for K and Wyj. The Gieren 
K data are a selected subset of the Persson et al. data dis- 
cussed next. The di sagreement in V may be a ttributed 
to instability width (fCa ldwcll k Coulson 198(|) and the 
placement of £ Car within that strip. 

6.3.2.3 - Persson et al. (2004) K-band: present the most 
extensive infrared photometry (CIT system) of LMC 
Cepheids. An infrared PLR is of interest because it is 
less sensitive to uncertainties in interstellar extinction 
and the intrinsic width of the relation is likely to be small. 
DC Their PLR can be written 

K = -3.261(/o 5 F - 1) + 12.775 (10) 

The standard error of the slope is 0.04. This slope is not 
significantly different from the one we determine from our 
galactic stars (-3.37±0.09, TabledU). From TableHSJwe 
find a K-band Galactic Cepheid zero-point -5.67±0.03 for 
a PLR with slope -3.26. This direct comparison yields a 
K-band LMC distance modulus of 18.45 ± 0.04, uncor- 
rected for metallicity effects. 

RP Using the method of reduced parallaxes and the 
Persson et al. slope, our parallaxes yield M(K) = - 
3.261(logP -1) -5.678 (± 0.033). We estimate that the 
uncertainty in the Persson et al. LMC zero-point (in the 
shorter period range where our Galactic stars are), is 
about 0.03 mag. Combined with our K-band PLR and 
its uncertainty the Persson et al. relationship yields an 
LMC modulus of 18.45 ±0.04 without metallicity correc- 
tion. We return to the metallicity issue in Section l"6. 4. II 

6.3.2.4 - OGLE (1999) W v/ : has produced the largest 
amount of LMC Cepheid photometry. In Figure [H 
we plot an apparent Wyj PLR for 581 Cepheids in 
the LMC. These data were carefully preened, selecting 
only Cepheids with normal light curves and amplitudes 
(|Ngeow et al. I [2001 . They were kindly communicated 
by G. Tammann, and provide the highly precise slope 
and zero point listed in Table [T5l 

DC Direct comparison of the Wyj and M W ^ VI ^ zero- 
points from Figures [Hand [5] yields an LMC distance mod- 
ulus 18.49 ±0.03 with no metallicity corrections. Con- 
straining the slope to the OGLE value results in the LMC 
distance modulus 18.51 ±0.04 listed in Table [TBI 
RP These data, when fit to the relation 

W(VI) = V- 2.45(V - I) = b(logP - 1) + A (11) 

yield slope, b = -3.29±0.01, which agrees within the un- 
certainties with the slope derived from our Galactic stars. 
For these we obtained b = -3.37 ± 0.12 and - 3.30 ± 0.12 
in two slightly different solutions (see Section 6.1 and 
Table [12]). Thus we find no evidence for a difference in 
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slope for the two galaxies. Nevertheless, because such 
a difference has been suggested (Sandage et al. 2004) it 
is desirable to compare Galactic and LMC Cepheids in 
the same period range. In the case of the Galactic stars 
this omits I Car. We therefore solve equation (11) for 
Cepheids with log P in the range 0.5 to 1.1. The OGLE 
LMC data then yield b = - 3.31 ± 0.03 and A = 12.64 ± 
0.01, not different from the values for the whole sample. 
Adopting a slope b = -3.31 we find a = -5.85 ± 0.04 in 
the equation 

M w(VI) =-3.3l(logP-l)+a (12) 

from our Cepheids. This leads to an LMC modulus of 
A-a = 18.50 ± 0.04, uncorrected for metallicity effects in 

Wy/. 

To summarize this section our PLR can be used to ob- 
tain LMC distance moduli by comparing (DC) our ab- 
solute zero-points with apparent magnitude zero-points 
from OGLE, Perrson et al. , and Gieren et al. . Con- 
straining the slopes of the PLR in Figure [5] to those de- 
termined from the Gieren and OGLE PLR in Figures [5] 
and [9] yields the zero points, a(DC), found in Table [T5l 
Comparing the zero-points for all four PLR we find LMC 
distance moduli 18.45 ±0.04 for the K band, 18.51±0.04 
for Wy/, 18.52 ±0.06 for the V band, and 18.49 ±0.03 
for the OGLE Wy/. RP yields the zero-points a(RP) in 
TableH5]and LMC distance moduli 18.50±0.04 for Wy/ 
and 18.42 ± 0.04 for K. These moduli remain uncorrected 
for possible metallicity effects. 

6.4. Extragalactic Distances: Applying Our PLR 

In this section wc apply our PLR to the LMC and NGC 
4258, comparing our derived distances with those from 
other investigators. In the case of the LMC we briefly 
describe our adopted metallicity corrections. 

6.4.1. Metallicity Corrections and LMC Distance Modulus 

Note that none of the LMC distance moduli de- 
rived a bove (Table [151) hav e metallicity corrections ap- 
plied. IMacri et all pOOei demonstrate that a metal- 
licity correction is necessary by comparing metal-rich 
Cepheids with metal-poor Cepheids in NGC 4258. 
With a previously m easured [O/H] metallicity gradient 
(|Zaritskv et al. Ill994[ ) Macri et al. find a Cepheid metal- 
licity correction in Wy/, 7 = —0.29 ± 0.09 r ± 0.05 s mag- 
nitude for 1 dex in metallicity, where r and s subscripts 
signify random and systematic. This value is similar 
to an earlier Wy/ metallicity correction (Kennicutt ct 
al. 1998) derived from Cepheids in M101 (-0.24± 0.16). 
Other less direct determinations (based for instance on 
RGB tip distances and Baade-Wesselink type luminosi- 
ties) are summarized by Macri et al. and are in agree- 
ment with these figures. Taking the weighted mean of 
the Kennicutt and Macri values and using the difference 
in metallicity of LMC and Galactic Cepheids (-0.36 dex 
from means of the data in Groenewegen et al. 2004 tables 
3 and 4) we find a metallicity correction of -0.10± 0.03 
magnitude with the Galactic Cepheids being brighter. 
The LMC distance moduli in Table [15] from the Pers- 
son et al. K data in the previous subsection suggest that 
the metallicity correction may be smaller for K than for 

Wy/. 

Returning to the issue of the true distance modulus to 
the LMC, our lowest error estimate is derived from the 



the OGLE photometry (Section 6.3.2.4, OGLE: m-M = 
18.50± 0.04). Combined with the estimated metallicity 
correction (-0.10± 0.03 m agnitude) we obtain an LMC 
modulus of 18.40 ± 0.05. iBenedict et al. I (|2002aD lists 
84 determinations complete through 2001 which can be 
compared with our new modulus. One recent determi- 
nation is noteworthy fo r its lack of dependence on any 
metallicity corrections. iFitzpatrick et all (|2003T l derive 
18.42± 0.04, from eclipsing binaries, a modulus in excel- 
lent agreement with our new value. 

6.4.2. NGC 4258 Distance Modulus 

Using HS T photometry of NGC 4258 Cepheids, 

IMacri et al. I (|2006l ) have determined a distance modu- 
lus relative to the LMC. They find that the difference 
in distance moduli of the LMC and NGC4258 is 10.88± 
0.06 mag. NGC 4258 has an independently measured 
distance obtained by comparing circ umnuclear maser 
prope r motions and radial velocities (jHerrnstein et al~l 
Il999| ). Macri et al. surveyed two fields in NGC 4258, one 
near the nucleus, the other in the outer regions of the 
galaxy. Macri (private communication) has applied our 
Wy/ PLR (Figure El) directly to N=85 inner field, so- 
lar metallicity Cepheids and finds m-M = 29.21 ± 0.02. 
The maser distance modulus is m-M = 29.29 ± 0.15. 
Our adopted LMC modulus, m-M= 18.40 ± 0.05 and 
the Macri differential modulus (LMC-NGC 4258) leads 
to a modulus of 29.28 ± 0.08 for NGC 4258, a value in 
even closer agreement with the maser-based distance. 

7. SUMMARY 

1. HST astrometry has now yielded absolute trigono- 
metric parallaxes for 10 Cepheid variables with an 
average av = 0.2 mas, or a^/n — 8%. These 
parallaxes, along with precision photometry culled 
from the literature, Lutz-Kelker-Hanson bias cor- 
rections, and reddening corrections derived from 
both the literature and our ancillary spectropho- 
tometry, provide absolute magnitudes with which 
to form Period-Luminosity relations. At logP = 
1, our zero-point errors are now at or below 0.03 
magnitudes in all bandpasses. 

2. Including perturbation orbits in our astrometry for 
W Sgr and FF Aql results in Cepheid orbit size 
and perturbation inclination. Assuming masses for 
the secondaries consistent with their known spec- 
tral type yields relatively low precision masses for 
these two Cepheids. We find for W Sgr and FF Aql, 
respectively, M = 6.5 ± 2M and M = 4.5 ± 1M . 
The major contributor to the mass uncertainty re- 
mains the parallax for FF Aql and the secondary 
spectral type for W Sgr. 

3. Comparing our parallax of T Vul with a parallax 
determined through the surface brightness tech- 
nique for that Cepheid, we find agreement at the 
la level. Alternately, if we demand that the surface 
brightness parallax be the same as our HST paral- 
lax, we determine a quasi-geometrical value of the 
radial velocity p-factor, p = 1.19 ± 0.16. Our PLR 
in the V magnitude agrees within la in slope and 
zero point with the Galactic PLR relation based on 
the Bayesian surface brightness PLR of Barnes et 
al. 
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4. Comparing our Wy/ zero-points with those 
adopted by the Freedman and Sandage Ho projects, 
we find better overall zero-point agreement with 
Sandage. However, the PLR adopted by Freedman 
et al. agrees with ours at longer periods. Given that 
most of the Cepheids observed in external galax- 
ies are long-period, there may be little effect on 
the Freedman et al. H value. Adopting our PLR 
would increase the Sandage et al. H value. 

5. Comparing our V, K, and W V i PLR with LMC 
PLR we find slope agreement for K and Wvi within 
the errors. The disagreement in V may be at- 
tributed to instability width and the placement of 
I Car within that strip. Comparing (both di- 
rect comparisons and via reduced parallaxes) zero- 
points yields a Wyi LMC distance modulus, (m- 
M) = 18.50 ±0.04, without any metallicity correc- 
tion. Adopting a metallicity correction -0.10± 0.03 
magnitude between Galactic and LMC Cepheids 
(with Galactic being brighter), we find a true LMC 
distance modulus (m-M) = 18.40 ± 0.05. 

6. Applying our PLR directly to Cepheids in NGC 
4258 provides a distance modulus, m-M = 29.21 
± 0.02, in good agreement with the maser distance 
modulus, m-M = 29.29 ± 0.15. From a metallicity- 
corrected LMC distance modulus and the Macri 
et al. difference in distance moduli we obtain m-M 
= 29.28 ± 0.08. 



LKH bias correction. We thank Lucas Macri for valu- 
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Log of Observations, Pulsational Phase, and Estimated B-V 
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a B-V estimated from phased light curve. 
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TABLE 2 

Target Cepheid Apparent Properties 



ID 


P(days) 


log P 


(V) 


(I) a 


(K) b 


(B-V) 


E(B-V) 


A v 


Ak 


£ Car 


35.551341 


1.5509 


3.732 


2.557 


1.071 


1.299 


0.17 


0.52 


0.06 


C Gem 


10.15073 


1.0065 


3.911 


3.085 


2.097 


0.798 


0.018 


0.06 


0.01 


Dor 


9.842425 


0.9931 


3.751 


2.943 


1.944 


0.807 


0.044 


0.25 


0.03 


W Sgr 


7.594904 


0.8805 


4.667 


3.862 


2.796 


0.746 


0.111 


0.37 


0.04 


XSgr 


7.012877 


0.8459 


4.556 


3.661 


2.557 


0.739 


0.197 


0.58 


0.07 


Y Sgr 


5.77338 


0.7614 


5.743 


4.814 


3.582 


0.856 


0.205 


0.67 


0.07 


8 Cep 


5.36627 


0.7297 


3.960 


3.204 


2.310 


0.657 


0.092 


0.23 


0.03 


FF Aql 


4.470916 


0.6504 


5.372 


4.510 


3.465 


0.756 


0.224 


0.64 


0.08 


T Vul 


4.435462 


0.6469 


5.752 


5.052 


4.187 


0.635 


0.064 


0.34 


0.02 


RT Aur 


3.72819 


0.5715 


5.464 


4.778 


3.925 


0.595 


0.051 


0.20 


0.02 



a Cousins I b CIT K 



TABLE 3 

t Car Reference Stars: Visible and Near-IR Photometry 



ID 


FGS ID 


V 


B-V 


U-B 


V-I a 


K a 


J-K a 


V-K a 


4273957 b 


2 


14.32 


0.71 


0.30 


0.89 


12.52 


0.30 


1.80 


4273905 b 


4 


13.53 


0.95 


0.63 


1.13 


11.20 


0.62 


2.33 


2M C 


5 


13.23 


1.18 


1.08 


1.33 


10.29 


0.78 


2.95 


4066585 b 


8 


10.77 


1.58 


1.95 


1.87 


6.56 


1.10 


4.22 


4066439 b 


9 


13.49 


0.57 


0.08 


0.72 








4066556 b 


10 


13.01 


0.60 


-0.04 


0.80 


11.48 


0.34 


1.53 


a Cousins I, Bessell/Brett JHK b 


ID from UCAC2 catalog c ID= 


09454541-6230004 from 2MASS catalog 







TABLE 4 

£ Car Field Av from Reference Star Spectrophotometry 



ID 


Av(B-V) 


Av(V-I) 


Av(V-K) 


Av(J-K) 


Av(U-B) 


<Av> 


SpT 


2 


0.4 


0.6 


0.4 


-0.3 


0.6 


0.33±0.19 


GOV 


4 


0.0 


0.5 


0.2 


0.2 


-0.2 


0.14 0.13 


G8III 


5 


0.5 


0.8 


0.7 


0.9 


0.7 


0.72 0.07 


K0III 


8 


0.6 


1.1 


1.0 


1.3 


1.2 


1.02 0.14 


K4III 


9 


0.5 


0.6 


0.0 


0.0 


0.2 


0.45 0.12 


F3V 


10 


0.6 


0.8 


0.5 


0.5 


-0.2 


0.46 0.18 


F4V 


(Av) 


0.45 


0.73 


0.55 


0.52 


0.38 


0.52 




± 


0.10 


0.09 


0.15 


0.31 


0.24 


0.06 





TABLE 5 

I Car Astrometric Reference Star Spectrophotometry Parallaxes 



ID 


V 


Sp. T. 


My 


Av 


m-M 


Tabs (mas) 


2 


14.32 


GOV 


4.4 


0.33 


9.89±0.5 


1.2±0.3 


4 


13.53 


G8III 


0.9 


0.14 


12.63 0.7 


0.3 0.1 


5 


13.23 


K0III 


0.7 


0.72 


12.48 0.5 


0.4 0.1 


8 


10.77 


K4III 


-0.1 


1.02 


10.72 0.5 


1.1 0.3 


9 


13.49 


F3V 


3.2 


0.45 


10.24 0.5 


1.1 0.3 


10 


13.01 


F4V 


3.3 


0.46 


9.67 0.5 


1.4 0.3 
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TABLE 6 

I Car and Reference Star Relative Positions a 



FGS ID V i b r] b 

£ Car 3.72 51.0107±0.0005 29.4377±0.0003 

2 14.32 20.4338 0.0006 133.6164 0.0005 

4 13.53 -61.6015 0.0007 64.6358 0.0007 

5 13.23 262.9969 0.0008 57.4027 0.0006 
8 10.77 199.6315 0.0006 36.9602 0.0004 
9 C 13.49 0.0000 0.0006 0.0000 0.0005 
10 13.01 161.3411 0.0006 -31.8256 0.0005 



a epoch 2004.431 b £ and rf are relative positions in arcseconds c RA = 09 45 07.44 , Dec = -62 30 57.9, J2000, epoch 2004.431 



TABLE 7 

i Car and Reference Star Relative Proper Motions 



ID Mx a vy a 



1 Car -0.0126±0.0003 0.0085±0.0004 

2 -0.0098 0.0010 0.0094 0.0011 

4 -0.0009 0.0012 0.0096 0.0013 

5 -0.0073 0.0017 0.0068 0.0016 

8 -0.0056 0.0012 0.0040 0.0011 

9 -0.0106 0.0009 0.0064 0.0008 

10 -0.0066 0.0010 0.0055 0.0011 



a fi x and fiy are relative motions in arcsec yr 



TABLE 8 

t Car and Reference Star Parallaxes and Transverse Velocities 



ID 


M a,b 


Tabs b 


Vt c 




(mas yr -1 ) 


(mas) 


(km s — 1 ) 


t Car d 


15.2 


2.01±0.20 


36±4 


2 


13.5 


1.33 0.14 


48 6 


4 


9.6 


0.32 0.04 


144 75 


5 


10.0 


0.45 0.05 


105 14 


8 


6.8 


1.19 0.10 


27 4 


9 


12.4 


0.73 0.22 


80 25 


10 


8.6 


1.45 0.11 


28 3 



a fi = (fi x + fiy) 1 / 2 from Table [Tp Final from modeling with equations 2 - 

5 C Vt = 4.74 X fi/Tr a bs d Modeled with equations 2-5, constraining D=-B and 
E=A 
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ASTROMETRIC REFERENCE STAR FINAL PARALLAXES 



ID 1 


V 


Sp. T. 


7r abs (mas 


LC-2 


14.29 


GOV 


1.3±0.1 


LC-4 


13.53 


G8III 


0.3 0.1 


LC-5 


13.18 


K0III 


0.4 0.1 


LC-8 


10.62 


K4III 


1.2 0.1 


LC-9 


13.44 


F3V 


0.7 0.2 


LC-10 


12.97 


F4V 


1.5 0.1 


ZG-2 


13.78 


G8III 


0.3 0.1 


ZG-3 


11.47 


F3.5V 


2.2 0.1 


ZG-5 


12.36 


F6V 


1.9 0.1 


ZG-8 


7.55 


G3V 


27.2 0.2 


ZG-10 


14.25 


F5V 


0.8 0.1 


ZG-11 


12.56 


KOIII 


0.5 0.1 


BD-2 


15.84 


AOIII 


0.1 0.1 


BD-3 


13.26 


F5V 


1.3 0.1 


BD-4 


15.79 


G3V 


0.6 0.1 


BD-5 


14.70 


G9V 


1.7 0.1 


BD-6 


15.28 


KOIII 


0.1 0.1 


BD-7 


15.29 


G5V 


0.9 0.1 


BD-8 


16.42 


KOV? 


0.7 0.1 


WS-4 


11.25 


F1V 


2.5 0.1 


WS-5 


13.25 


KOIII 


0.9 0.2 


WS-7 


12.8 


KOIII 


0.6 0.1 


WS-9 


14.17 


F8V 


1.5 0.1 


WS-10 


13.7 


MOIII 


0.3 0.1 


WS-11 


14.1 


F2III 


0.3 0.1 


XS-2 


14.00 


KOIII 


0.5 0.1 


XS-3 


13.10 


B7V 


0.5 0.1 


XS-4 


13.62 


Aim 


0.5 0.1 


XS-5 


12.56 


KOIII 


0.9 0.1 


XS-6 


13.04 


F5V 


1.9 0.1 


XS-7 


12.56 


F3V 


1.9 0.1 


XS-8 


13.98 


A1V 


0.7 0.1 


YS-2 


10.37 


A5V 


2.2 0.3 


YS-3 


12.41 


A5V 


1.0 0.1 


YS-4 


13.36 


KOIV 


1.6 0.2 


YS-7 


11.18 


FOV 


2.2 0.2 


YS-9 


14.92 


K7V 


4.9 0.5 


YS-10 


12.83 


G9III 


0.5 0.1 


FF-2 


14.17 


K2III 


0.3 0.1 


FF-3 


14.16 


K3V 


3.6 0.2 


FF-4 


13.68 


K3V 


4.0 0.2 


FF-5 


14.93 


G7V 


1.6 0.1 


FF-6 


15.1 


F2V 


0.6 0.1 


FF-7 


15.29 


K2III 


0.2 0.1 


TV-2 


13.79 


KOIII 


0.4 0.1 


TV-3 


13.31 


G3V 


2.1 0.2 


TV-4 


14.29 


K1IV 


0.7 0.1 


TV-5 


13.26 


GOV 


1.5 0.2 


TV-6 


11.69 


K1.5III 


0.8 0.1 


TV-7 


14.48 


KOIV 


0.6 0.1 


TV-8 


12.60 


K3III 


0.6 0.1 


RT-4 


13.87 


K2V 


2.7 0.2 


RT-5 


13.26 


KOIII 


0.4 0.1 


RT-6 


11.37 


G2V 


3.0 0.3 


RT-7 


11.47 


F3V 


2.4 0.2 


RT-8 


13.90 


F3V 


0.9 0.1 


RT-9 


14.93 


G5III 


0.2 0.1 



1 LC-2 = i Car, Ref-2; ZG = C Gem; BD = beta Dor; WS = W Sgr; XS = 
X Sgr; YS = Y Sgr; FF = FF Aql; TV = T Vul; RT = RT Aur 
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TABLE 10 
Binary Cepheid Orbits and Masses 



Parameter 


W ogr 


r r Aql 


d(mas) 


2.67 ± 0.2 


3.36 ± 0.4 


Pfdavs) 


1582 ± 3 


1434 ± 1 


P ( vears i 


4.33 ± 0.01 


3.93 ± 0.01 


To 


2004.16 ± 0.01 


2003.29 ± 0.04 


e 


0.41 ± 0.02 


0.09 ± 0.01 


i 


7?0 ± 0?8 


33°± 5° 


n 


68?4 ± 4?0 


61.3°± 9° 




328? ± 1?3 


327° ± 4° 


Secondary Sp. T. 


A5V-F5V 1 


F1V 2 


Secondary Mass 


2.0-1.4 


1.6 


a (mas) 


12.9 ± 0.3 


12.8 ± 0.9 


a (AU) 


5.67 ± 0.13 


4.54 ± 0.14 


f 


0.207 ± 0.017 


0.263 ± 0.031 


Cepheid Mass (Mq) 


6.5 ± 2 


4.5 ± 1 



1 Range from N. Evans (private communication) 2 lEvans "eTaTI p99"0h 



TABLE 11 

Cepheid Parallaxes, Proper Motions, and Absolute Magnitudes 



Parameter 






Cepheid 








t Car 


C Gem 


/3 Dor 


W Sgr 


X Sgr 


Duration (y) 


1.95 


1.50 


1.50 


1.71 


1.49 


Ref stars (#) 


6 


6 


5 


6 


7 


Ref (V) 


13.00 


12.03 


14.95 


13.04 


13.28 


Ref (B-V) 


0.92 


0.69 


0.77 


1.35 


0.98 


w abs (mas) 


2.01 ±0.20 


2.78 ±0.18 


3.14 ±0.16 


2.28 ±0.20 


3.00 ±0.18 


H (mas y" 1 ) 


15.2±0.5 


6.2±0.5 


12.7±0.8 


6.6±0.4 


10.0±1.2 


PA. (°) 


304±2 


272±5 


10.4±0.6 


134±8 


193±3 


Av 


0.52±0.06 


0.06±0.03 


0.25±0.05 


0.37±0.03 


0.58±0.1 


LKH Corr 


-0.08 


-0.03 


-0.02 


-0.06 


-0.03 


(m-M)o 


8.56 


7.81 


7.50 


8.31 


7.64 


M v 


-5.35±0.22 


-4.03±0.15 


-4.05±0.11 


-3.97±0.20 


-3.68±0.17 


M/ 


-6.31±0.22 


-4.80±0.15 


-4.74±0.11 


-4.62±0.20 


-4.64±0.17 


M K 


-7.55±0.21 


-5.73±0.14 


-5.62±0.11 


-5.51±0.19 


-5.15±0.13 


M W(VI) 


-7.71±0.21 


-5.92±0.14 


-5.76±0.11 


-5.58±0.19 


-5.28±0.13 




Y Sgr 


5 Cep 


FF Aql 


T Vul 


RT Aur 


Duration (y) 


1.51 


2.44 


1.77 


1.67 


1.49 


Ref stars (#) 


6 


5 


6 


6 


5 


Ref (V) 


12.51 


12.06 


14.48 


13.39 


13.02 


Ref <B-V) 


0.98 


1.30 


1.16 


1.12 


0.80 


w abs (mas) 


2.13 ±0.29 


3.66 ±0.15 


2.81 ± 0.18 


1.90 ±0.23 


2.40 ±0.19 


fi (mas y" 1 ) 


7.0±0.8 


17.4±0.7 


7.9±0.8 


7.1±0.3 


15.0±0.4 


PA. (°) 


204±5 


-73±3 


144±11 


141±6 


179±3 


Av 


0.67±0.04 


0.23±0.03 


0.64±0.06 


0.34±0.06 


0.20±0.08 


LKH Corr 


-0.15 


-0.01 


-0.03 


-0.12 


-0.05 


(m-M)o 


8.51 


7.19 


7.79 


8.73 


8.15 


M v 


-3.42±0.30 


-3.47±0.11 


-3.05±0.15 


-3.24±0.28 


-2.90±0.18 


M 7 


-4.12±0.30 


-4.14±0.11 


-3.65±0.27 


-3.79±0.28 


-3.47±0.18 


M K 


-5.00±0.30 


-4.91±0.09 


-4.39±0.14 


-4.57±0.24 


-4.25±0.17 




-5.04±0.30 


-5.09±0.10 


-4.53±0.13 


-4.69±0.27 


-4.37±0.17 
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TABLE 12 

Galactic Cepheid PLR Zero-points (a) and Slopes (b) 



Source 1 


B07 


B07f 


F01 


S04 


B03 


a 

V 


-4.05±0.02 


-4.03±0.03 


-4.22+0.02 


-4.00± 0.10 


-4.16± 0.22 


I 


-4.78 0.03 




-4.90 0.01 


-4.78 0.01 




K 


-5.71 0.03 


-5.67 0.02 








W VI 


-5.85 0.03 


-5.83 0.03 


-5.90 0.01 


-5.96 0.04 




b 












V 


-2.43 0.12 


-2.44 0.11 


-2.76 0.03 


-3.09 0.09 


-2.69 0.17 


I 


-2.81 0.11 




-2.96 0.02 


-3.35 0.08 




K 


-3.32 0.12 


-3.35 0.08 








Wyj 


-3.37 0.12 


-3.30 0.12 


-3.26 0.01 


-3.75 0.12 





1 B07 = this paper; B07f = no LKH, this pape r; F01 = IFreedman et al. I j200lT> ; S04 = 
ISandage eTTaD H2004I) . B03 = IBarnes etajTI II2OO30 . All PLR are parameterized as M = a 
+b*(logP-l). 



TABLE 13 

Reduced Parallax (RP) Zero-point Tests of Freedman and Sandage 



Source 


ZP 


F (Eq. 8) 


5.823 ± 0.036 


Freedman 


5.979 


Difference 


+0.156 


S (Eq. 9) 


5.964 ± 0.042 


Sandage 


5.959 


Difference 


-0.005 



TABLE 14 
Cepheid Parallaxes and Postdictions 



ID 


HST 7r ai)s (mas) 


S06 a 


Diff 


F01 b 


Diff 


I Car 


2.01±0.2 


1.74 


0.27 


1.88 


0.13 


C Gem 


2.78 0.18 


2.73 


0.05 


2.64 


0.14 


P Dor 


3.14 0.16 


2.94 


0.2 


2.84 


0.3 


W Sgr 


2.28 0.2 


2.34 


-0.06 


2.19 


0.09 


XSgr 


3 0.18 


2.9 


0.1 


2.69 


0.31 


YSgr 


2.13 0.29 


2.02 


0.11 


1.84 


0.29 


8 Cep 


3.66 0.15 


3.96 


-0.3 


3.6 


0.06 


FF Aql 


2.81 0.18 


2.68 


0.13 


2.39 


0.42 


T Vul 


1.9 0.23 


1.88 


0.02 


1.68 


0.22 


RT Aur 


2.4 0.19 


2.4 





2.11 


0.29 


unwt'd mean Diff(10) 






+0.052 




+0.225 


std err 






±0.052 




±0.039 


std dev 






0.155 




0.140 



a Parallax predicted from Equation 9 and Sandage zero-point b Parallax predicted from Equation 8 and Freedman zero-point 



TABLE 15 

LMC PLR Zero-points (a). Slopes (b), and Distance Moduli 



Source 1 


Band 


a 


b 


a(DC) 2 


a(RP) 3 


DC(m-M) 4 


RP(m-M) 4 


G05 


V 


14.42±0.05 


-2.78±0.09 


-4.10±0.04 




18.52 ±0.06 






K 


12.78 0.02 


-3.26 0.03 


-5.70 0.03 




18.48 0.04 






Wyj 


12.65 0.02 


-3.37 0.03 


-5.86 0.03 




18.51 0.04 




Per04 


K 


12.78 0.02 


-3.26 0.03 


-5.70 0.04 


-5.68 0.03 


18.48 0.04 


18.45 0.04 


OGLE 


W V j 


12.65 0.01 


-3.29 0.01 


-5.84 0.03 


-5.85 0.04 


18.49 0.03 


18.50 0.04 



1 G05 = IGieren et al. I | |2005T) ; Per04 = "P ersson et al. I ||2004T) ; OGLE = private communication G. Tammann. All PLR are parameterized 
as M = a +b(logP-l). 2 Zero-points obtained by fitting the data plotted in Figure[5]but with slopes constrained to those from G05, Per04, 
and OGLE. 3 Zero-points obtained via reduced parallaxes. 4 Distance moduli with no metallicity corrections applied. 
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Fig. l — 



I Car field with astrometric reference stars marked. Box is 15' across. 
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0.4 - 



0.6 - 



0.8 - 



1.0 - 




i.2 -q 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 

1 2 3 4 5 

V-K 

Fig. 2. — J-K vs V-K color-color diagram for I Car and reference stars. The dashed line is the locus of dwarf (luminosity class V) stars 
of various spectral types; the dot-dashed line is for giants (luminosity class III). The reddening vector indicates A^=1.0 for the plotted 
color systems. For this low-Galactic latitude field (Ay) = 0.52 ± 0.06 magnitude (TableQ. 
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J-K (2Mass) 

Fig. 3. — Reduced proper motion diagram for 430 stars in a 1/3 degree field centered on I Car . Star identifications are shown for I 
Car ('1') and our astrometric reference stars. % for these stars is calculated using our final proper motions (TableO. For a given spectral 
type giants and sub-giants have more negative Hj{ values and are redder than dwarfs in J-K. Reference stars ref-4, -5, and -8 are confirmed 
to be giant stars. The cross in the lower left corner indicates representative errors along each axis. 
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Fig. 4. — Histograms of x and y residuals obtained from modeling i Car and astrometric reference stars with equations 4 and 5, 
constraining D=-B and E=A. Distributions are fit with gaussians whose l-o dispersions arc noted in the plots. 
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Fig. 5. — Period-Luminosity relations for V, I, K, and Wesenheit Wvi, where Wvi = (V)-2.45((V)-(I)). Coefficients are for M 
= a + b * (logP — 1). Errors are 1<t. 



Galactic Cepheid Period-Luminosity Relations 



23 



> 




-4 ■ 



1.0 1.2 
logP 

Period-Luminosity relations for the present work compared to those adopted by Frccdman ct ai. (2001) (left) and 
ISandage et ai. I ((2004) (right). For the Freedman comparison we use 2.45 as our Wv; color coefficient. For the Sandage comparison we 
recalculate our Wvi using a 2.52 coefficient. While the Sandage et aJ. PLR agrees better with the bulk of our Cepheids, the Freedman 
et a!. PLR better matches ours for logP>l. Most Ho Key Project extragalactic Cepheids are long-period. Hence our new PLR has minimal 
impact on the Freedman et ai. Hq value. If adopted, our PLR would tend to increase the Sandage et ai. (2006) Hq value. 
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Fig. 7. — Period-Luminosity relation for My where we fit the seven Cepheids in our sample with logP< 1. This fit excludes i Car , f 
Gem, and /? Dor which has logP = 0.9931. Coefficients are for M = a + b(logP — 1.0). The slopes differ by ~ 2cr, comparing the logP < 1 
slope against the entire range (Figure Qj). Our parallaxes provide only weak evidence for a break at logP = 1. 
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Fig. 8. — Period-Luminosity relations for absorption-corrected V, K, and Wyj for 13 LMC Cepheids from Giercn ct al. (2005). The 
K-band magnitudes have been transformed to the CIT system. Coefficients are for M = a + b(logP — 1.0). The K-band and Wyi slopes 
are in better agreement with the Galactic Cepheid results than the V-band slope (Figure [5j- 
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Fig. 9. — Wyj Period-Luminosity relation for 581 LMC Cepheids from OGLE. Coefficients are for M = a + b(logP — 1.0) and Wyj 
<V)-2.45«V)-<I». 



